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ABSTRACT. In this work, we study superconvergence properties for some high-order 
orthogonal polynomial interpolations. The results are two-folds: When interpolating func- 
tion values, we identify those points where the first and second derivatives of the inter- 



polant converge faster; When interpolating the first derivative, we locate those points where 
the function value of the interpolant superconverges. For the earlier case, we use various 
Qh' Chebyshev polynomials; and for the later case, we also include the counterpart Legendre 



polynomials. 
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1. Introduction 

In numerical computation, we often observe that the convergent rate exceeds the best 
possible global rate at some special points. Those points are called superconvergent points, 
\ and the phenomenon is called superconvergence phenomenon, which is well understood 

for the h- version finite element method, see, e.g., [H HH H21 [261 EZl EH EH [39l SHI E2] 
and references therein. As comparison, the relevant study for the p-version finite element 
method and the spectral method is lacking. Only very special and simple cases have been 
^ \ discussed in the following works: [321 150} [5T]. 

CN | The study of the superconvergence phenomenon for the /i-version method has made 

great impact on scientific computing, especially on a posteriori error estimates and adaptive 
methods, which is well documented in the following books: [U El l27l [46] and works cited. It 
is the believe of this author that the scientific community would also be benefited from the 
study of superconvergence phenomenon of spectral collocation methods as well as related 
^-version and spectral methods. This work is the first step, where the superconvergence 
points of some orthogonal polynomial interpolation will be identified. 

The most celebrated advantage of spectral methods is the exponential (or geometric) rate 
of convergence for sufficiently smooth, essentially analytic functions. However, most error 
bounds in the literature are in the form of iV _fc ||zi||jjfc+i, where ./V is the polynomial degree 
(or trigonometric function degree in the Fourier spectral case). We can see this in almost 
all books on spectral methods such as [3 [3 (H El [TO1 [HI [HI [HI H3 [231 IM1 [251 E21 SH H2] 
Ideally, we expect to establish the convergence rate p~ N for some p > 1 or e~ aN for some 
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a > 0. There have been some limited discussion of this type of error bounds in the past, 
e.g., in these two books |19^ 147] . see also \20\ [371 145| . In the framework of p- and hp- finite 
element methods, results of the exponential convergent rate can be found in books \40\ |4"3"] 
and references therein, see also, \21\ I22j. 

Actually, exponential rate of convergence for polynomial approximation of analytic func- 
tions can be traced back 100 years. The following result is due to Bernstein [6], see also 
[30|, Theorem 7] : / is analytic on [— 1 , 1] if and only if 

sup hm tfE n {f) = p~\ E n (f)= inf \\f - g^, 

n->oo g&Vn 

where V n is the polynomial space of degree no more than n, and p > 1 is the sum of the 
half-axes of the maximum elliptic disc D p bounded by the ellipse E p with foci ±1 that / 
can be analytically extended to. This fact serves as a starting point of the first part of 
our analysis in this paper. We shall focus our attention to the approximation properties of 
polynomial interpolants of analytic functions on [—1,1]. We identify those points in [—1, 1], 
where the derivatives are superconvergent, in the sense that the convergent rate gains at 
least one factor iV _1 . Three popular Chebyshev polynomial interpolants are discussed here: 
Chebyshev, Chebyshev-Lobatto, and Chebyshev-Radau. 

In the second part of analysis, we consider polynomial interpolants of the first deriva- 
tive of a one-order higher polynomial, and identify those points where function values are 
superconvergent. In this part, both the Chebyshev and Legendre polynomials are studied. 

Throughout of the paper, we use the standard notation for orthogonal polynomials: Tj, 
the Chebyshev polynomial of the first kind; Uj, the Chebyshev polynomial of the second 
kind; and Lj, the Legendre polynomial. 

2. Interpolation by the Chebyshev Polynomials 

In this Section, we discuss interpolations by the Chebyshev polynomials of the first kind, 
Chebyshev-Lobatto polynomials, and Chebyshev-Radau (right and left) polynomials. 

2.1. Statement of the Results 

We state the superconvergence results in this subsection. The meaning of the supercon- 
vergence will be made clear in the next subsection. 
2.1.1. Chebyshev interpolant. Interpolating at the zeros of Tjv+i, 

2k 4- 1 

r^+i(x fc ) = cos(iv + i)e fc = o, x k = cose k , e k = ^—TT, k = o,i,...,N. (2.1) 

Proposition 2.1. For the Chebyshev interpolant, the first derivative superconverges at 
zeros of Un, which are 

^ = COS iVTI' fc = l,2,...,JV; (2.2) 
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and the second derivative superconverges at cos 9 k with 9 k satisfies the following equation: 

(N + 1) cos(iV + 1)0 sin 9 = sin(iV + 1)9 cos 9. (2.3) 

This set of cos9k are close to interior zeros (not including k = 0, N in (|2.ip ) of Tjy + i for 
large N. 

2.1.2. Chebyshev-Lobatto interpolant. Interpolating at the zeros of Tjv+i — Tjy-i, 

(T N+1 -T N ^)(x k ) = 2smN9 k sm9 k = 0, x k = cos^, k = 0, 1, . . . , N. (2.4) 
We also call this set of points the 2nd type Chebyshev points. 

Proposition 2.2. For the Chebyshev-Lobatto interpolant, the first derivative supercon- 
verges at cos 9k with 9k satisfies the following equation: 

N cos N9 sin 9 + sin N9 cos 9 = 0. (2.5) 

2k — 1 

This set of #&s are close to zeros of cos N9 for large N, i.e., 9k ~ ~2N~ 7T w ^ en ^ k * s awa y 
from and tt. 

As for super convergent points of the second derivative, 9 k s satisfy the following equa- 
tion: 

(iV 2 - 1) sin N9 sin 9 = 2N cos N9 cos 9. (2-6) 

This set of 9kS are close to interior zeros (not including and tt) of sin N9 for large N, i.e., 
9k ~ — , when 9k is away from and tt. 

Remark 2.1. We see that when interpolating at the zeros of T/v+i(x), the derivative of the 
interpolant superconverges at the zeros of Un(x), and the second derivative super convergent 
points go back "almost" to the zeros of T/v+i(x) except the two ends; when interpolating 
at the zeros of (T/v + i — T/v-i)(x) = 7jv(1 — x 2 )Un-i(x), the derivative of the interpolant 
superconverges "almost" at the zeros of Tn(x), and the second derivative superconvergent 
points go back "almost" to the zeros of Un-i(x). 

2.1.3. Chebyshev-Radau (right) interpolant. Interpolating at the zeros of T/v+i — T/v, 

1 9i 2kir 

(T N+1 -T N )(x k ) = 2sm(N + -)9 k sm^- = 0, x k = cos 2jV+r k = 0,1, . . . ,N. (2.7) 

Proposition 2.3. For the right Chebyshev-Radau interpolant, the first derivative super- 
converges at cos 9k with 9k satisfies the following equation: 

19 19 

(2N + 1) cos(iV + -)9 sin - + sin(iV + -)9 cos - = 0. (2.8) 

2k — 1 

This set of 9kS are close to zeros of cos(iV + ^)9 for large N, i.e., 9k ~ ^ — ^ tt, when 9 k 
is away from 0. 
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As for superconvergent points of the second derivative, 9k satisfies the following equation: 

19 19 

(2N 2 +2N + 1) sm(N + -)9 sin - = (2N + 1) cos(iV + -)9 cos -. (2.9) 

This set of O^s are close to interior zeros (not including 0) of sin(A r + \)9 for large N, i.e., 
2/c7r 

9u , when 9h is away from 0. 

27V + 1 

2.1.4. Chebyshev-Radau (left) interpolant. Interpolating at the zeros of T/v+i + T/v, 

(T N+1 + T N )(x k ) = 2cos(iV + ^)6» fe cos y = 0, x fc = cos + ^ = 0, 1, . . . , iV. 

(2.10) 

Proposition 2.4. For the left Chebyshev-Radau interpolant, the first derivative supercon- 
verges at cos 9k with 9k satisfies the following equation: 

16 19 

(2N + 1) sm(N + -)9 cos - + cos(iV + -)0 sin - = 0. (2.11) 

1 2k 
This set of #&s are close to zeros of sin(iV + -)9 for large N, i.e., 0& « — ^ ir, when 0^ 

is away from ir. 

As for superconvergent points of the second derivative, 6k satisfies the following equation: 

16 16 

(2N 2 + 2N + 1) cos{N + -)6 cos - = (2JV + 1) sin(iV + -)9 sin -. (2.12) 

This set of 9kS are close to interior zeros (not including tt) of cos(A r + \)9 for large N, i.e., 

(9z. ~ — — 7r, when 9u is away from 7r. 

2iV + l 

Remark 2.2. Propositions 2.3 and 2.4 say that when interpolating at the right (left) 
Radau points, derivative of the interpolant super converges "almost" at the left (right) Radau 
points except x = —1 {x = 1), and the second derivative superconvergent points "almost" 
go back to the interior right (left) Radau points. 

Remark 2.3. In this section, we provide derivative super convergence points as roots 
of some polynomial equations, and the approximated values of those points. Numerical 
data in the last section demonstrate that these approximation are usually good enough for 
practical purpose. Therefore, they can be used as the intitial guesses, e.g., the Newton 
iteration, should the more accurate values are desired. On the other hand, there are quick 
and simple ways to find roots of sums of orthogonal polynomials by expressing the problem 
as an eigenvalue problem wuith the matrix being upper Hessenberg, see |15| . 

2.2. Analysis 

Let u be analytic on / = [—1,1]. According to Bernstein [BJ, u can be analytically 
extended to D p , which is enclosed by an ellipse E p with ±1 as fosi, p > 1 as the sum of its 
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semimajor and semiminor: 

E p : z = ^(pe ie + p- l e~ ie ), 0<6<2ir. 

We consider polynomial G T'jv who interpolates it at N + 1 points — 1 < xq < x\ < • • • < 
xtv < 1. The error equation is, according to [131 P-68], expressed as 

1 f u N+1 (x) u(z) ^ ^_„ttJV 



it(x) - Ujv(x) = — T-ir dz, u N+l (x) = cU j=0 (x - xj). (2.13) 

2m J Ep ujn+i[z) z-x J 

By direct differentiation of (|2.13p . c.f., [37], one obtains the error equation for the derivative 
and the second derivative 

«.(,) _ „ M ,2.J E ^ + f^ + 1^>) ( , 15) 

We need to estimate wjv+i(x), a;^ v+1 (x), and u>^ +1 (x) on [—1, 1] to establish the error 
bounds. However, by each differentiation, we lose at least one power of N. 

Key observation. Let us examine the error equation (|2.14p . At the N special points 
u'n+i( x ) = Oj we have only the second term, which is usually smaller than the first term in 
magnitude by a factor N or TV 2 as we will see later. Similarly, at the N — 1 special points 
u n+i( x ) = 0; we have only the second and third terms left in the error equation (j2. 15|) . 
Again, we may gain a factor N in the error bounds. 

We consider the four sets of interpolation points in the previous subsection. 

The exponential decay of the error is provided by the value of lon+i(z) on the ellipse 
E p in the denominators of ()2.13j) - fj2.15p . In all four sets of interpolation points, ujn + ±(z) 
involves the Chebyshev polynomials of the first kind. We have the following characteristic 
expressions. The proof is elementary by the definition of and therefore is omitted. 

Lemma 2.1. For z £ E p , we have 



\T N+1 (z)\ = p 2N+2 + p -2N-2 + 2cos2 (N + 1)6- (2.16) 

\T N+l (z) - T N - X (z)\ = i VP 2 + P~ 2 ~ 2 cos 26 ^ p 2N + p~ 2N - 2 cos 2N6; (2.17) 

\T N+l (z) ± T N {z)\ = -y/l + p~ 2 ± 2p~ l cos 6^p 2N+l + p~ 2N ± 2pcos(2AT + 1)6. (2.18) 
We also need to bound the derivatives of ujn+i(x) for x G [—1,1]. 
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Lemma 2.2. 



max \T' N+1 (x)\ 
xel-i,i] 


= (N + l) 2 , 


(2.19) 


max \T' N+l (x) - Tjy^x)] 
e[-i,l] 


= AN, 


(2.20) 


max |T^ +1 (x)±7^(x)| 

x6[— 1,1] 


= 2N 2 + 2N + 1; 


(2.21) 


max T^ +1 (x) 

£•€[-1,1] 


= ^N(N + 1)\N + 2), 


(2.22) 


max |T^ +1 (x) - 3^_i(a;)| 


= lAT(2iV 2 + l), 


(2.23) 


max |T^ +1 (x) ±T^(x)| 
xe[— 1,1] 


= ^JV(iV + l)(iV 2 + iV + l). 


(2.24) 



The proof is also elementary and is provided in the appendix. To establish the error 
bounds, we define some constants, which include C p (u) = max|u(z)|, D p , the shortest 

z^Ep 

distance from E p to [—1, 1], and L p , the arch length of the ellipse E p . We have 

D P = ^(P + P' 1 )-1, L p < ttVp 2 + P" 2 . (2.25) 

The latter is the Euler's estimate which overestimates the perimeter by less than 12 percent. 
Theorem 2.1. Assume that u is analytic on [—1, 1] and can be extended analytically to 
the region bounded by an ellipse E p . Let 6 Vn\~ 1j 1] be the interpolant of u at N + 1 
zeros of Tn+i(x) on [—1,1]. Then 

4Si " wWI - C 'X' ^ • (2 ' 26) 



C p (u)L p / AT(iV + l) 2 (iV + 2) 2(iV + l) 2 _2_\ 1 
I 3D P + D 2 + u» J 

(2.28) 

Furthermore, at those special points where T' N+1 (x) = 0, we have 



max \u"{x)-u%{x)\ < - ^ I 1J + ^ t ^ + 

xG[-l,l] 7T 



i // \ / C p (u)L p 1 J7T 

rng^ |n fa) - u N {tj)\ < ^ _ , t^cos^^, j = 1,2 

and at those special points where Tj^ +1 (x) = 0, we have 



i • • • ? ' j 



AT, 
(2.29) 



where Tj = cos #j with 9j satisfies ()2.3[) . 



6 



Proof: Since un £ 7-V[— 1, 1] interpolates u at zeros of Tjv+i(x), we have wtv+i = TXr+i 
in (|2.13p - (|2.15p . By the identity (j2. 16|) . we can derive the lower bound 

l^iWI^f^ 1 -^- 1 ). 

Substituting this into f|2. 13j) and using max \Tn + i(x)\ = 1, we derive, 

se[-i,i] 

\u{x)-u N {x)\<- I ^\d\z\ 1 Vx€[-l,l]. (2.31) 

Note that |it(z)| < C p (u) and |z — re | ~ 1 < D" 1 for z £ E p and x E [—1, 1], and we obtain 
O from (IOTP . 

Using (|2.19p in (|2. 14[) and following the same procedure as above, we derive (|2.27p . 
Similarly, we establish (1238]) by applying (1239]) and (12321) in (pT[5D . 

At the special points when T^ +1 (x) = (iV + 1)Un(x) = 0, the first term on the right- 
hand side of (|2.14p is gone, we then obtain (|2.29p following the same argument as we derive 

(23511 ■ 

At the special points when T^ +1 (x) = (N + l)U' N (x) = 0, the first term on the right- 
hand side of ()2.15p is gone, we then obtain (I2.30P following the same argument as we derive 
(|2,27p . Since Un(x) has N simple roots in (—1, 1), it is guaranteed that U' N {x) has N — 1 
simple roots in between and they can be expressed as xt = cos 6 k with 6k satisfying (I2.3j) . 



Theorem 2.2. Under the same assumption as in Theorem 2.1, let un G Vn[—1, 1] be the 
interpolant of u at N + 1 zeros of Tn + \(x) — T^-i{x) on [—1, 1]. Then 



C p (u)L p (p-p- 1 )- 1 
x6[-l,l] vrZ) p p w - p- N 



max \u(x) - u N (x)\ < n — ^ rxr, ( 2 - 32 ) 



max - uV*)| < (V + ±) ^ P , (2.33) 

xe[-i,i] kD p V A*/ P ~P 

max |</'f*) - «« fx) I < Cp{u)Lp (" 4iV(2Ar2 + 1) + 8iV + M f 2 34) 

Furthermore, at those special points where T^ +1 (x) — T^_ 1 (x) = 0, we have 

max Wit,) - u' N ( h )\ < C i^ pi ^^ p r N \ ( 2 -35) 

where tj = cos6j with 6j satisfies ()2.5|) ; and at those special points where Tjy +1 (x) — 
2~iv-i( x ) = ^' we have 



max 

i<j<N-i 



I /// \ It t \\<T 2C p( U ) L P f A AT, 1 \ jP~P ^ 1 /O 

l« fa) - ujvfa)| < ^ y m+ D~ p ) P N - P ~ N ' (2,36) 
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where Tj = cos 9j with 9j satisfies (|2.6p . 

Proof: Since un G T^Arf— 1,1] interpolates u at zeros of (Tv+i — T/v-i)(a;), we have 
ajjv+i = Tjv+i — in (|2.13p - f|2. 15|) . By the identity (|2.17p . we can derive the lower 

bound 

\(T N+1 - T N ^){z)\ > \{p - p^){p N - p~ N ). 
Using (|2.20p and ()2.23|) . the rest is then similar to the proof of Theorem 2.1. 

Theorem 2.3. Under the same assumption as in Theorem 2.1, let itjv G Vn[—1, 1] be the 
interpolant of it at iV + 1 zeros of Tn+i(x) ± T^{x) on [—1, 1]. Then 

■IKi - ""Wl * °it' ^-7" ' (2 - 37) 



IKu (l) - u " <x)l - { — so, — + m + m) P » - p-" • (2 - 39) 



xe[-i,i] 

Furthermore, at those special points where T' N+1 (x) ± T' N (x) = 0, we have 

j^w^-Mm < C it^ [ )-- P ly ^ (2 - 40) 

where tj = cos8j with Oj satisfies ()2.11|) in case T/v+i + T/v, and (j2.8j) for T/v+i — T/v; and 
at those special points where T^ +1 (x) ± T'^{x) = 0, we have 

-1W 



max 

1<j'<JV-1 



VXr,) - <(r,)\ < « (4AT + i-) fi^., (2.4!) 



where Tj = cos 0j with 0j satisfies (|2.12p in case Tv + i + T/v, and (|2.9p for Tjv+i — Tv. 

Proof: Since it/v G P^f— 1, 1] interpolates it at zeros of (Tv+i ±T\r)(x), we have wat+i = 
T/v+i ± Tv in ()2.13p - (|2.15p . By the identity (|2.18p . we can derive the lower bound 

|(Tv +1 ± T N )(z)\ > i(l - p^){ P N - p~ N ). 
Using (|2.2ip and (I2.24|) . the rest is then similar to the proof of Theorem 2.1. 

We see from the above that superconvergent points involve extremals of ujn+i{x) for 
x G [— 1, 1]. It would be interesting to see the distribution and magnitudes of those points 
for different cases. Since those information is clear for the case un+i{x) = T/v +1 (x), we 
consider the other three cases. 

Theorem 2.4. The envelope for the extremals of T/v+i — T/v-i on [—1, 1] forms an ellipse 

y 2 

x 2 H — 7T = 1 (see Figure [T]). 
2^ 



8 



Proof: We need to demonstrate that the envelope of (x, Tjv+i (x) — T/v-i(x)), or 



(cos 9, cos(N + 1)0 - cos(N - 1)0) or (cos 0, -2 sin NO sin 9) 

2j — l 

is an ellipse. Note that siniV0 has extremals at 9j = ir and points (cos 9j, ±2sin0j) 

are on the indicated ellipse. 

Theorem 2.5. The envelope for the extremals of T/v+i ± TW on [—1, 1] form a parabola 
2(l±x) = y 2 (see Figure |5J). 

Proof: We need to demonstrate that the envelope of (x,T/v+i(x) ± Tjy(x)), or 

(cos 0, cos(N + 1)6* ± cos N9), 

or 

19 19 

(cos0,2cos(iV + -)0cos-) and (cos 9, -2sin(iV + -)0sin -) 

are parabola. Note that 

1 1 
cos(iV + -)6l and sin(iV+-)0 

2 2 



have extremals at 



2jV 2j + l 

and p,- = 7r, 



J 2iV + l J 27V + 1 

respectively, and points 

9 ■ 9 ■ 

(cos 9j, ±2 cos -^) and (cos 6>j, ±2 sin -^-) 

are on the indicated parabola, respectively. Indeed, 

9 9 
2(1 + cos 9) = (2 cos-) 2 and 2(1 - cos 9) = (2 sin -) 2 . 



3. Derivative Interpolation 

In this part, we consider different interpolants to the first derivative of a smooth function, 
and identify superconvergent points for the function value approximation. To be more 
precise, we construct polynomial un € Vn such that 

u N (-l) = u(-l), u' N (x k ) = u'(x k ), k = l,2,...,N, -1 < x x < ■ ■ ■ < x N < 1, (3.1) 

and locate the point yk, where (u — u^)(yk) is superconvergent. It is worthy to point 
out that this supeconvergence knowledge can be utilized in spectral collocation method for 
solving ODEs. We shall demonstrate this point later in our numerical examples. 

To fix the idea, we consider only the case u € Vn+1, since superconvergent property 
may be narrowed down to the capability of a polynomial space to approximate polynomials 
of one order higher [4] . 
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In addition to Chebyshev, Chebyshev-Lobatto, and Chebyshev-Radau interpolants, we 
also consider Gauss, Gauss-Lobatto, and Gauss-Radau interpolants. 
We begin with two technical lemmas. 

Lemma 3.1. 

N(L N + L N _t)(x) = (x + 1)(L N - L N ^)'(x). (3.2) 

N(L N - L N - X ){x) = (x- 1)(L N + L N -i)'{x). (3.3) 

Proof: We only prove (|3.2p . Multiplying both sides of (|3.2p by x — 1 and using the 
identity 

(x 2 - l)L' N {x) = {N + l)L N+1 (x) - (N + l)xL N (x), 

we have 

N(x-l)(L N + L N -i)(x) 
= (N + l)L N+ i(x) - (AT + l)xL N {x) - NL N {x) + NxL N _i(x). 

Canceling and collecting the same terms on two sides, we obtain 

(2N + l)xL N (x) - iVLjv-i(x) = (N + l)L N+1 (x), 

which is the three-term recurrence relation. 

Remark 3.1. We see that zeros of the derivatives of the right (left) Legendre-Radau 
polynomials are zeros of the left (right) Legendre-Radau polynomials. 

Lemma 3.2. 

(T N + T N -!)(x) = (x + 1)(C/jv-i - U N . 2 )(x). (3.4) 

(T N - T N -!)(x) = (x- 1)(C/jv-i + U N . 2 )(x). (3.5) 
Proof: By the definition, 

1 9 

(T N + T N _i)(x) = cos NO + cos(iV - 1)6 = 2cos(iV - -)6>cos-, 

TT , sin NO sm(N-l)9 cos(N - ±)0 

[U N -i - U N - 2 ){x) = — — — = - e . 

sin 6 sin 6 C os | 

Therefore, 

(T N + T N ^)(x) = 2cos 2 -(C/jv_i -U N - 2 ){x) 

= (1 + cos0)([/ 7V -i - «7jv_ 2 )(a;) = (1 + x)(C/jv-i - C/jv_ 2 )(x). 

This established (|3.4p . The proof of (|3.5p is similar and hence is omitted. 
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Remark 3.2. There is an interesting similarity between (|3.4p and (|3.2p . as well as (|3.5p 
and (|3.3p . Note that T' N {x) = NUn-i(x), and therefore, zeros of the derivatives of the right 
(left) Chebyshev-Radau polynomials are 'almost' zeros of the left (right) Chebyshev-Radau 
polynomials. 

Theorem 3.1. Let u G Vn+i an d G "Pat satisfy (|3.ip with the collocation points x^s 
being the roots of L^. Then we have (up to a constant), 

f x 1 x 2 -1 

(u - u N )(x) = J ^ L N (t)dt = m - - (L N+1 - L N -x)(x) = N , N + js L N (x). (3.6) 

Proof: We see that (u — un)' G Vn- By the definition of the interpolation points, we have, 
up to a constant, (u — un)'(x) = Ljy(x). Using the initial condition itjv(— 1) = u(— 1), the 
conclusion follows by integration. 

Remark 3.3. Theorem 3.1 says that when interpolating derivative at the N Gauss points, 
the function value approximation is superconvergent at the N — 1 interior Lobatto points, 
i.e., roots of L' N . 

Theorem 3.2. Let u G Vn+i and un G Vn satisfy (|3.ip with the collocation points x^s 
being the roots of Ln — Ln-2- Then we have (up to a constant), 

" = NfNTT) jV (A r + 1)(2jV _3) ^- 1 " Ln -^' ^ 

Proof: From (u — un)' G Vn, we have (up to a constant), (u — un)'{x) = Ln{x) — Ln-2(x). 
By the identity 

[(1 - x 2 )L' n (x)}' + n(n + l)L n (x) = 0, 

we have 

u>(x) u> (x) _ [(l-^_ 2 (x)r [(1 - x 2 )L' N {x)\ 
u \ x ) u n\ x ) — (A r_ 2)(A r_i) N{N + 1) ' 

Integrating with the initial condition and using the identity 

(2n + = (-L n +i - ^n-l)'(^), 

we derive 

(l-x 2 )L5v- 2 (^) (l-x 2 )L' N (x) 



u(x) — un(x) 



(JV-2)(JV-1) iV(iV + l) 



(x 2_ 1)(2Ar _ 1)Ljv i(x) 4i y_ 2 

JV(JV + 1) (N + 1)N(N- l){N-2)) [ J % - 2i J ' 

which is the right-hand side of (j3.7|) by the identity 

1 : (L n +l(x) - Ln-iix)) = 1 (x 2 - l)L' n {x). 

in + 1 n(n + 1) 
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Remark 3.4- We see that the magnitude of the first term on the right-hand side of 
(|3.7p is larger than that of the second term by a factor about N. Therefore, the function 
value approximation reaches its best at roots of Ljv-i- In other words, Theorem 3.2 says 
that when interpolating derivative at the N Lobatto points, function value approximation 
is super convergent at the N — 1 Gauss points. 

Theorem 3.3. Let u € Vn+i and 6 Vn satisfy f)3. 1 [) . and let the collocation points XkS 
be the roots of =p Ljv-i- Then we have (up to a constant), 

iV 2 N 
(u - u N )(x) = N 2 _i ( Ln ^ l n-i)(x){x T 1) - N 2 _i ( Ln T ^JV-i)(a;)(a; ± !)• (3.8) 

Proof: By (|3.3|) . we have (up to a constant), 

(u - u N )'(x) = N(L N =F L N -i){x) = {L N ± L N ^ 1 )'(x)(x =F 1). 

Integrating both sides with the initial condition and using f)3.2|) - (|3.3|) . we have 

(u - u N )(x) = J i L N± L N ^)'(t){t T l)dt 

= {L N ± L N - X )(x)(x Tl)-^{L N T L n ^){x){x ± 1) + J { l nT L N - X ){t)dt. 

Multiplying both sides by N and moving the last term on the right-hand side to the left-hand 
side, we obtain (|3.8p . 



Remark 3. 5. We see that the magnitude of the first term on the right-hand side of (|3.8h is 
larger than that of the second term by a factor N. Therefore, function value approximation 
reaches its best at the roots of Ln ± Lat_i for (I3.8p . In other words, Theorem 3.3 says 
that when interpolating derivative at the N left (right) Radau points, the function value 
approximation is superconvergent at the N right (left) Radau points. 

Now we turn to the Chebyshev polynomials. We need the following two identities: 

(VT^T^x))' + —L=T n {x) = 0, (3.9) 
yl-r 

-(Tn-^x) - T n+l {x)) = (x 2 - l)U n -i(x), (3.10) 

Theorem 3.4. Let u 6 Vn+i an d £ Vn satisfy (I3.ip . and let the collocation points XkS 
be the N roots of T^. Then we have (up to a constant), 

i \( \ N ( 2 -im ( \ xT N {x) + {-l) N 

(u-u N )(x) = -p— - {x -l)U N -i{x) N2 _ . (3.11) 

Proof: We have (up to a constant) (n — un)'{x) = Tn(x). By (|3.9p . we write 



(« - u N y(x) = -^L^(^T^r N (x)y. 
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Integrating and using the initial condition, we derive 



(u - u N )(x) = J ^ T N (t)dt = -J N 2 (Vl-t 2 T' N {t))'dt 
= -^T' N (x) - ^ JjT' N (t)dt 

x 2 -i TT xT N ( x ) + (-1)^ i 

N 2 

Moving the last term to the left and multiplying the resultant by — ^ , we obtain (|3.11jl . 

Remark 3.6. We see that the magnitude of the first term on the right-hand side of 
(|3.1ip is larger than that of the second term by a factor TV. Therefore, the function value 
approximation reaches its best at roots of C/jv-i- I n other words, Theorem 3.4 says that 
when interpolating derivative at the N Chebyshev points of the first kind, function value 
approximation is superconvergent at the N — 1 Chybyshev points of the second kind. 

Theorem 3.5. Let u G Vn+i and un £ Vn satisfy (13. ip . and let the collocation points XkS 
be ±1 plus the N — 2 roots of Un-2- Then we have (up to a constant), 

Nil-x 2 )^ . , N(l-x 2 ) TT . . xT N (x) + (-l) N 

iv 2 (xr^ 2 (x) + (-i)") jv-i r 

+ 2{N-2) 2 {N 2 -l) + (iV2-l)(iV-2)2 ^W*- I 3 - 12 ) 

Proof: We have (up to a constant), (u — un)'(x) = (x 2 — 1)Un-2(x). By (13.10P and 
(|3.9p . we write 

(u - u N )'(x) = ^(T N _ 2 (x) - T N (x)) 

( 1 , n w , w 1 



(^(V^^rO*))' - ^l_^ (v /T3^_ 2 (*))' 



2 

Integrating and using the initial condition, we derive 



2 



X 



+^^|z^ + £ IW-W* + ^ £pw-w - ?V(*))«. 
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N 2 

Moving the last term to the left-hand side and multiplying the resultant by —5 , we 

N z — 1 

derive (|3.12p when replacing Ujy-i — Un-3 by 2Tjv-i- 

Remark 3. 7. We see that the magnitude of the first term on the right-hand side of (|3,12|) 
is larger than that of the rest terms by a factor about N. Therefore, (u — un)(x) reaches its 
best at roots of Tjv-i- In other words, Theorem 3.5 says that when interpolating derivative 
at ±1 plus the N — 2 Chebyshev points of the second kind, function value approximation 
is superconvergent at the N — 1 Chebyshev points of the first kind. 

Theorem 3.6. Let u 6 Vn+i an d un £ Vn satisfy (13. ip . and let the collocation points XkS 
be the roots of T/v ± Tjv-i- Then we have (up to a constant), 

w x N{x±l) /rT1 m w . N(x 2 -l) , . xT N (x) + (-l) N 

(u-u N )(x) = -^^(T NT T N _ 1 )(x)± {N2 \ ){N ) _ 1) U N - 2 (x)T N{ J 2 _\ ] 

— {N 2 -i)(N-iy — ± (iv 2 _ i)(jv-i)2 T ^) dt - ( 3 - 13 ) 

Proof: The proof is similar to that of Theorem 3.5 and hence is omitted. 

Remark 3.8. We see that the magnitude of the first term on the right-hand side of 
(|3.13p is larger than that of the rest terms by a factor about N. Therefore, the function 
value approximation reaches its best at the roots of T/v =F Ijv-i for (|3.13|) . In other words, 
Theorems 3.6 says that when interpolating derivative at the iV left (right) Chebyshev-Radau 
points, function value approximation is superconvergent at the N right (left) Chebyshev- 
Radau points. 

4. Numerical Tests 

In this section, we perform numerical tests on two typical analytic functions. 

Example 1. f(x) = (1 + 25x 2 )~ 1 . This is the well known Runge's example [35] • Deriva- 
tive errors of its interpolants at the Chebyshev points, the 2nd type Chebyshev points (or 
Chebyshev-Lobatto points), and the right-Chebyshev-Radau points are depicted in Figure 
[3l Figure [H and Figure [5l respectively. Derivative superconvergent points are marked by 
'*'. We see that the errors at the superconvergent points are significantly smaller (by a 
magnitude) than the maximum error just as Theorems 2.1 - 2.3 predicted. 

Next, we solve the following initial value problem 

1 1 

u (x) = T7, u(—l) = — 

v ; l + 25x 2 ' v ; 26 

by collocating at the Chebyshev points as in (4.1). From Figure [6l we see that when 
interpolating the derivative at the Chebyshev points, the function value of the interpolant 
converges much faster at the interior 2nd type Chebyshev points marked by '*', which is 
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Figure 1: Profile of the Chebshev-Lobatto Polynomial 

in consistent with Theorem 4.1. We also see that the function value approximation has 
maximum error at the derivative interpolation points marked by 'o'. 

Analytic function f(z) = (1 + 25X 2 )" 1 has two single poles at ±i/5 and it is straightfor- 
ward to calculate p = (\/5 2 + 1 + l)/5 ~ 1.2198. Therefore, we expect a slow convergence 
as we have observed form Figure [3] - Figure El 

Example 2. f(x) = (2 — an analytic function which has a simple pole at x = 2 

and p = 2 + \/2 2 — 1 ~ 3.7321. Therefore, we expect much faster convergence compared 
with Example 1. It is indeed the case. We plot counterparts graphs in Figure 0— Figure 
[TO] and we see that high accuracy is achieved with relative very low polynomial degree n. 
We observe similar super convergence phenomena as in Example 1. 

Conclusion Remarks. The results in Section 2 can be extended to all Legendre based 
polynomials as in Section 3. The reason for using the Chebyshev based polynomials is that 
they have simple (trigonometric functions) expressions. 

Extension of the result to more general Jacobi type polynomials is feasible. However, 
the analysis would be much more involved. 
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